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DEFINITION  OF  SYMBOLS  (Concluded) 
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cylinder 
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respectively 


KyXr  specific  twist  of  ring 

A  =  qm2 

^x>  Poisson's  ratios  associated  with  bending  in  the  x  and  y 

directions,  respectively 

H i  Poisson's  ratios  associated  with  extension  in  the  x  and  y 

x’  y  directions,  respectively 

cr  longitudinal  stress  in  corrugation 


vii 


THE  GENERAL  INSTABILITY  OF  RING-STIFFENED 
CORRUGATED  CYLINDERS  UNDER  AXIAL  COMPRESS  ION 

SUMMARY 


■r* 

method  is  presented  to  determine  the  general  instability  load  of  a  ring- 
stiffened  corrugated  cylinder  under  axial  compression.  This  method  is  developed 
using  linear  small  deflection  theory.  The  stiffness  properties  of  the  rings  are 
uniformly  distributed  along  the  cylinder  and  the  eccentricity  of  the  rings  with  re¬ 
spect  to  the  corrugation  centerline  is  taken  into  account^ 


Analytical  and  experimental  results  are  compared.  In  this  comparison 
good  agreement  is  obtained  for  cylinders  loaded  in  pure  compression.  For  the 
cylinders  subjected  to  bending  or  a  combination  of  bending  and  compression  the 
analytical  calculations  are  conservative. 


A  computer  program  for  the  application  of  this  method  has  been  developed. 
The  program  and  an  explanation  of  its  notation  are  included  in  this  report. 


INTRODUCTION 


Axially  loaded  ring-stiffened  corrugated  cylinders  are  susceptible  to  two 
types  of  instability;  panel  instability  and  general  instability.  Panel  instability  ijsf 
defined  as  buckling  of  the  cylinder  between  rings.  In  general  instability,  the  rings 
as  well  as  the  cylinder  undergo  buckling  displacements.  Panel  instability  is  dis¬ 
cussed  in  many  textbooks  and  manuals;  however,  the  same  is  not  true  for  general 
instability. 

In  comparing  the  stiffness  characteristics  of  a  corrugated  cylinder  to  those 
of  a  monocoque  cylinder  of  the  same  weight  it  must  be  noted  that  the  corrugated 
shell  has  a  much  greater  longitudinal  bending  stiffness.  This  advantage  is  offset 
to  some  degree  by  the  low  circumferential  extensional  stiffness  of  the  corrugated 
cylinder.  Because  of  this  lack  of  circumferential  extensional  stiffness  each 
corrugation  will  act  as  an  individual  column  unless  the  shell  is  reinforced  by  rings. 
Low  extensional  stiffness  may  be  desirable,  however,  in  areas  of  high  thermal 
gradients. 

In  the  past,  ring-stiffened  cylinders  have  often  been  designed  to  prevent 
general  instability  by  employing  semi-empirical  methods  to  size  the  rings.  One 
of  the  most  commonly  used  of  such  methods  is  given  by  Shanley  [  1] .  Until  re¬ 
cently  it  was  believed  that  the  ring  moment  of  inertia  obtained  from  his  formula 


was  sufficient  to  safeguard  against  general  instability.  Tests  performed  under 
the  direction  of  the  Marshall  Space  Flight  Center  as  part  of  the  Saturn  V  develop¬ 
ment  program  have  shown,  however,  that  Shanley's  criterion  may  be  very  uncon¬ 
servative  for  the  design  of  ring- stiffened  corrugated  cylinders.  The  results  of 
these  tests  are  given  in  Tables  I  and  II. 


Since  Shanley's  work  was  published,  considerable  research  has  been 
performed  on  the  stability  of  stiffened  cylinders.  A  major  contribution  was  made 
by  Van  der  Neut  [  2]  ,  when  he  considered  the  effect  of  ring  and  stringer  eccentric¬ 
ities.  Additional  contributions  were  made  by  Hedgepeth  and  Hall  [3] ,  Card  [4], 
and  Baruch  and  Singer  [5] .  Often,  however,  research  in  this  field  has  produced 
a  method  which  is  either  too  academic  or  too  complex  to  be  used  by  the  stress 
anal\ 


The  purpose  of  this  report  is  to  present  a  reliable  and  relatively  simple 
meansijQf  predicting  the  general  instability  load  of  corrugated  cylinders  under  axial 
compression.  This  method  considers  the  eccentricity  of  the  rings  with  respect 
to  the  corrugation  centerline.  It  also  incorporates  all  the  stiffnesses  attributable 
to  the  rings  and  the  shear  and  longitudinal  stiffness  properties  of  the  corrugation. 
The  circumferential  extensional  and  bending  stiffnesses  and  the  torsional  stiffness 
of  the  corrugation  are  small  (generally  less  than  1  percent  of  the  longitudinal 
stiffnesses)  and  are  therefore  neglectecLj| 


ns 
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GENERAL  THEORY^ 
Assumptions 


||The  method  of  analysis  presented  in  this  paper  is  based  on  the  following 
assumptions: 


1.  Linear  small  deflection  theory  applies.  This  is  justified  because 
the  high  longitudinal  bending  stiffness  of  a  corrugated  cylinder  makes  it  less 
susceptible  to  initial  imperfections  and  other  monocoque  effects. 


2.  The  longitudinal  wavelength  of  the  buckled  skin  is  sufficiently  large 
to  permit  a  "smeared"  ring  approach;  i.  e.  ,  all  the  ring  stiffness  parameters 
may  be  uniformly  distributed  along  the  cylinder. 

3.  The  corrugated  cylinder  can  be  treated  as  an  equivalent  orthotropic 

cylinder,  the  radius  of  which  is  equal  to  the  mean  radius  of  the  corrugation.  / 

-» . — 4 
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The  circumferential  extensional  and  bending  stiffnesses  of  the  cor¬ 
rugation  as  well  as  its  torsional  stiffness  can  be  neglected  since  they  are  small 
as  compared  to  the  longitudinal  stiffness  properties  of  the  shell. 

5.  Buckling  displacements  are  sinusoidal. 


6.  Prebuckling  deformations  are  neglected. 


Plasticity  effects  and  local  failures  are  not  considered,  jj _ _ 


/  o 


Displacements  and  Boundary  Conditions 


The  cylinder  is  in  equilibrium  Under  the  applied  load  just  prior  to  buckling 
and  deformations  due  to  buckling  are  measured  from  this  position.  In  accordance 
with  assumption  5,  displacements  may  be  written  in  the  form 


u  =  U  cos  — —  cos  n0 

i-l 

T  T  .  1X1 7TX  •  ,  j 

v  =  V  sin  — —  sin  n0  (1) 

i-l 

.  m7TX 

w=  W  sin — ~ —  cos  n0 
L 

This  corresponds  to  the  following  simply  supported  boundary  conditions  at 
x  =  0,  L 

w  =  0  N  =  0 

x 

v  =  0  M  =  0 

x 

Thus,  at  the  ends  of  the  cylinder  motion  radially  and  tangentially  is  prevented, 
while  longitudinal  motion  is  allowed;  i.e.  ,  u  ^  0.  These  boundary  conditions  are 
appropriate  for  cylinders  bounded  by  deep  supporting  rings,  which  are  rigid  in 
their  own  planes  but  may  readily  bend  or  warp  out  of  their  planes. 


Elastic  Relations 


For  an  orthotropic  shell  the  relations  for  the  stress  resultants  and  couples 
in  terms  of  the  strains  and  curvatures  are 
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M  =  D  K 
xy  xy  xy 


In  view  of  assumption  (4)  and  with  the  use  of  the  reciprocal  theorem, 
equations  (2)  lead  directly  to  the  following  expressions  for  the  corrugated  cylinder 


=  =  3U 

N  =  E  e  =  E  ~r— 

X  x  ox 

Nxy  R  dd  +  3x 

M  =  D  k  =  D  |-t 
xxx  x  ox 


N  =  M  =M  =0 
y  y  xy 


(3) 


The  rings  are  considered  eccentric  with  respect  to  the  centerline  of  the 
skin  and  displacements  of  the  ring  at  its  centroid  may  be  found  by  means  of  the 
transformations 


u  =  u  -  c 
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9w 

ax 


V 

r 


_r  _c  9w 
R  V  R  80 


(4) 


w  =  w 
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where  c  is  the  distance  from  the  ring  centroid  to  the  skin  centerline.  For  a  ring 
loaded  normal  to,  as  well  as  in  the  plane  of,  initial  curvature,  the  strain,  curva¬ 
tures  and  specific  twist  at  the  central  axis,  may  be  written 
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Assuming  the  shear  center  to  coincide  with  the  ring  centroid,  the  stress 
resultants  and  stress  couples  acting  on  an  element  of  the  ring  are  related  to  the 
strains  and  curvatures  by  the  equations 
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yr 
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Substituting  equations  (4)  and  (5)  in  the  above  expressions  gives 
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Equilibrium  Equations 


(6  concluded) 


Differential  equations  are  obtained  by  considering  separately  the  equi¬ 
librium  of  an  element  of  the  skin  and  a  corresponding  "smeared"  ring  element. 

Six  conditions  of  equilibrium  must  be  satisfied  by  the  forces  and  moments  acting 
on  each  of  the  elements.  An  element  of  the  skin  is  shown  in  Figure  1  a,  b,  c. 

The  first  of  theso  figures  shows  the  stress  resultants  and  the  forces  X,  Y,  and  Z 
acting  at  the  interface  of  the  skin  and  the  ring,  and  the  second  figure  gives  the 
stress  couples  and  the  moments  Tx  and  Ty  transferred  into  the  skin  element  from 
the  ring.  Components  of  the  buckling  force  resulting  from  the  deformation  of  the 
element  are  shown  in  Figure  ic.  The  component  in  the  x  direction  is  multiplied 
by  the  quantity  (1  +  ex)  to  include  the  straining  of  the  middle  surface  as  suggested 
by  Fliigge  (7).  The  components  in  the  y  and  z  directions  are  due  to  the  rotations 
9v/ 8x  and  8w/ dx  of  the  element  respectively. 

The  six  conditions  of  equilibrium  for  the  skin  lead  to  the  following  equations 
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(7) 


These  equations  may  be  combined  and  put  in  the  form 
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The  six  equations  of  equilibrium  for  an  element  of  the  ring  (Fig.  2)  may 
be  obtained  in  a  similar  manner.  They  are 
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The  last  three  of  the  above  equations  are  used  to  eliminate  Qyr  and  Nyxr* 
This  yields  the  following  three  equations 
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The  forces  X,  Y,  and  Z  and  the  moments  Tx  and  Ty  acting  at  the  inter¬ 
face  of  the  skin  and  ring  may  now  be  eliminated  by  subtracting  equation  (8)  from 
equation  (10).  Hence, 
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These  are  the  equilibrium  equations  in  terms  of  the  stress  resultants. 
Substituting  the  expressions  (3)  and  (6)  for  the  stress  resultants  into  equations 
(11)  gives  the  equilibrium  equations  in  terms  of  the  displacements. 
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After  introducing  the  expression  for  the  displacements  given  by  equations 
(1)  into  the  differential  equations  (12)  one  obtains  the  following  three  linear  equa¬ 
tions  . 


+  q  m2  W  =  0 


Determination  of  Buckling  Load 

The  set  of  homogeneous  equations  (13)  has  nontrivial  solutions  only  when 
the  determinant  of  its  matrix  is  zero.  Equations  (13)  may  be  written  in  matrix 
form 


“an  +  A 

a12 

a2i 

a22  +  A 

a3i 

a32 

a13 

U 

a23 

V 

a33  +  A 

w 

—  M 

(14) 
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The  determinant  of  the  characteristic  matrix  of  A  is  a  polynomial  of  the 
third  degree  in  X.  The  problem  of  determining  the  buckling  load  (q  =  X/n?  )  for 
known  values  of  m  and  n  has  therefore  been  reduced  to  that  of  finding  the  roots 
(eigenvalues)  of  the  characteristic  equation 

a  +  bA  +  cA2  +  A3  =  0  (15) 

where  a,  b,  and  c  are  known  functions  of  the  coefficients  ay.  Obviously  only 
real  and  positive  roots  of  equation  (15)  are  of  interest,  the  lowest  of  which  de¬ 
termines  the  buckling  load  for  the  mode  shape  under  consideration.  The  critical 
buckling  load  of  the  cylinder  may  be  found  by  calculating  q  for  a  range  of  values 
of  m  and  n,  and  plotting  a  family  of  curves  as  shown  in  Figure  3.  The  critical 
buckling  load  will  then  be  the  minimum  value  of  q  corresponding  to  integer  values 
of  m  and  n.  If  the  computer  program  (see  Appendix)  is  used,  the  minimum 
buckling  load  will  be  indicated  for  the  specified  range  of  m  and  n,  and  no  plotting 
is  required. 

Since  W  is  usually  in  the  order  of  n  times  larger  than  V  or  U  an  approxi¬ 
mate  solution  for  q  may  be  obtained  by  dropping  the  terms  containing  q  in  the 
first  two  of  equations  (13).  This  yields 


q 

m,  n 

where  |A|  is 


= _ lA] _ _ 

(ai2  -  ail  a22)  na2 

the  determinant  of  the  symmetrical  matrix 


an 

a12 

a13 

a12 

a22 

a23 

a13 

a23 

a33 

(16) 


If  an  iterative  procedure  is  used  to  find  the  roots  of  equation  (15) ,  the  above 
approximation  is  found  to  be  very  useful  as  a  starting  point. 


/comparison  WITH 


TEST  RESULTS)) 


(Axial  LoadJ 


A 


^=2 


As  a  part  of  the  Saturn  V  development  program  eleven  ring- stiffened 
corrugated  cylinders  were  tested  to  failure  under  axial  load.  Of  these,  five 
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/"failed  in  general  instability.  The  remaining  cylinders  failed  either  by  panel  in¬ 
stability  or  local  crippling  before  the  general  instability  load  was  reached.  V  A 
photograph  of  a  typical  general  instability  failure  is  shown  in  Figure  4.  "Figure  5 
shows  one  of  the  rings  of  this  cylinder  after  the  failure. nTable  I  gives  jthe  geom¬ 
etries  of  these  five  cylinders  along  with  their  actual  and  calculated  failure  loacls^- 
Thpf^ose  agreement  between  the  experimental  data  and  calculations  supports  the 
validity  of  the  method  of  analysis  given  in  this  report.  [/  — — 

r-— -p,  • 

|  Stable  i.  ring-stiffened  corrugated  cylinders  -  axial  loading  j 

|  _  liar w" 


(7 


CYLINDER 

N>) 

No.  2 

No.  3 

No.  4 

No.  5 

Aluminum  Alloy 

7075-T6 

7075-T6 

7075-T6  x 

7075-T6 

7075-T6 

Cylinder  Length  (in. ) 

33.0 

33.  0 

33.  0 

69.  6 

69.6 

(cm) 

83.  8 

83.8 

83.  8 

176.8 

176.8 

Radius  (in. ) 

24.7 

24.7 

24.7 

49.4 

49.4 

(cm) 

62.7 

62.7 

62.  7 

125.  0 

125.  0 

Corrugation  Pitch  (in. ) 

1.  43 

1.43 

1.43 

2.85 

2.85 

(cm) 

3.63 

3.63 

3.  63 

7.24 

7.24 

Corrugation  Thickness  (in.) 

0.  020 

0.  019 

0.025 

0.041 

0.041 

(cm) 

0.  05 

0.  048 

a  06 

0. 10 

0. 10 

Corrugation  Depth  (in. ) 

0.44 

0.  44 

0.  44 

0.87 

0.87 

(cm) 

1.12 

1.  12 

1.12 

2.21 

2.21' 

Type  of  Ring 

c* 

I 

X 

I 

X 

Ring  Spacing  (in. ) 

6.  38 

6.38 

6.38 

12.4 

12.4 

(cm) 

16.21 

16.21 

16.21 

31.5 

31.5 

Ring  Moment  of  Inertia  (in. 4) 

0.  0050 

0.  0104 

0.0104 

0.286 

0.286 

(cm4) 

0.  208 

0.433 

0.433 

11.9 

11.9 

Ring  Area  (in. 2) 

0.  040 

0. 121 

0.121 

0. 180 

0.180 

(cm2) 

0.  258 

0.78 

0.  78 

1. 16 

1.16 

Ring  Eccentricity  (in. ) 

-0.73 

-0.  53 

-0.  53 

-1.99 

-1.99 

(cm) 

-1.  85 

-1.35 

-1.35 

-5.05 

-5.  05 

Actual  Failure  Load  (Kips) 

131. 

174. 

224. 

659. 

648. 

(N) 

5.  83  x  10s 

7.  74  x  105 

9.  96  x  106 

2.93  x  106 

2.  88  x  106 

Calculated  Failure  Load  (Kips) 

118. 

198. 

233. 

654. 

654. 

(N) 

5.25  x  105 

8.81  x  105 

1.04  x  10  5 

2.91  x  105 

2.91  X  106 

Percent  Error  (%) 

11.0 

-12.1 

-3.9 

0.8 

-0.9 

*  For  this  ring  an  effective 
moment  of  inertia  and  ef¬ 
fective  area  were  calculated 
using  the  approach  given  in 
Reference  9. 
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The -instrumentation  on  these  cylinders  provided  some  interesting  informa¬ 
tion.  The  radial  deflection  gages  located  on  the  rings  showed  that  as  the  load 
increased  most  of  the  rings  deflected  into  the  theoretically  predicted  circumferen¬ 
tial  mode  shape.  Figure  6  shows  a  typical  example  of  this  phenomenon.  When  the 
radial  deflection  at  a  point  of  maximum  bending  on  the  ring  is  plotted  versus  axial 
load  the  curve  obtained  is  asymptotic  to  the  failure  load  as  is  shown  in  Figure  7. 

Another  interesting  phenomenon  was  shown  by  the  strain  gages  located  on 
the  corrugated  skin.  Two  gages  were  located  opposite  each  other  at  various  points 
on  the  skin,  see  Figure  8,  so  that  longitudinal  bending  of  the  corrugation  could  be 
observed.  As  the  axial  load  on  the  cylinder  increased,  the  strain  versus  load  plot 
of  the  two  gages  in  each  set  remained  linear  and  almost  coincident  until  just  below 
the  failure  load.  Then  at  various  locations  on  the  cylinder  the  strain  in  one  gage 
increased  nonlinearly  while  the  strain  in  the  other  leveled  off  and  then  decreased  as 
shown  in  Figure  9.  This  meant  that  the  corrugation  was  bending  appreciably,  and 
failure  of  the  cylinder  was  imminent.  The  actual  cylinder  failures  were  similar  to 
monocoque  failures  in  their  swiftness.  Conjecture  can  be  made  that  this  second 
phenomenon  might  be  used  to  determine  the  general  instability  failure  load  of  a 
ring-stiffened  corrugated  cylinder  without  actually  failing  the  specimen. 


J^Test  data  is  also  available  for  a  ring-stiffened  corrugated  cylinder  loaded 
in  pure  bending  and  for  a  cylinder  loaded  simultaneously  in  bending  and  axial  com¬ 
pression.  Table  II  gives  the  cylinder  geometries  and  loads \  j  Cylinders  loaded  in 
the  above  manner  are  usually  analyzed  by  calculating  the  maximum  compressive 
stress  due  to  bending  and  axial  load  and  then  assuming  this  stress  to  act  along  the 
entire  periphery  of  the  cylinder.  This,  however,  leads  to  a  conservative  prediction 
for  the  buckling  stress,  since  a  stiffened  cylinder  can  withstand  a  greater  maximum 
stress  in  bending  than  in  pure  compression.  This  increased  load  carrying  capability 
of  cylinders  in  bending  can  be  expressed  in  terms  of  a  bending  factor,  /3.  If  the 
bending  factor  j8  and  the  allowable  stress  in  pure  compression  cr  are  known  then  the 
maximum  allowable  stress  in  any  combination  of  bending  and  compression  may  be 
calculated  from  the  equation 


[(Axial  Load  and  Bending  J 


<J 

max 


a 


rd-^) 
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(TABLE  II.  RING-STIFFENED  CORRUGATED  CYLINDERS  - 
_  ’  AXIAL  AND  BENDING  LOADING 

fi  _ _  '  l  Wn  C  l  Ma  7 


CYLINDERS 

No.  6 

No.  7 

Aluminum  Alloy 

7075-T6 

7075-T6 

Cylinder  Length  (in. ) 

268.6 

33.0 

(cm) 

682.2 

83.8 

Radius  (in. ) 

197.6 

24.7 

(cm) 

501.9 

62.7 

Corrugation  Thickness  (in. ) 

0.  145 

0.020 

(cm) 

0.368 

0.051 

Corrugation  Pitch  (in. ) 

11.40 

1.43 

(cm) 

28.95 

3.63 

Corrugation  Depth  (in. ) 

3.48 

0.44 

(cm) 

8.  84 

1.12 

Ty  pe  of  Ring 

I 

C* 

Ring  Spacing  (in. ) 

51.0 

6.38 

(cm) 

129.  54 

16.20 

Ring  In- Plane  Moment  of  Inertia  (in.4) 

34.4 

0.  0050 

(cm4) 

0.  1431 

0.208 

Ring  Out-of- Plane  Moment  of  Inertia  (in. 4) 

0 

0 

(cm4) 

0 

0 

Ring  Torsional  Stiffness/ GR  (in.4) 

0 

0 

(cm4) 

0 

0 

Ring  Area  (in.2) 

2.48 

0.040 

(cm2) 

16.  00 

a  258 

Ring  Eccentricity  (in. ) 

-6.24 

-0.73 

(cm) 

-15.  84 

-1.85 

Actual  Bending  Load  (Kip-ft. ) 

27,  900 

161.0 

(N-m) 

3.79  x  107 

2.  18  x  105 

Actual  Axial  Load  (Kips) 

6,  930 

0 

(N) 

30. 8  x  106 

0 

Actual  Maximum  Stress  (lb.  / in. 2) 

43,  530 

37,  960 

(N/m2) 

3.00  x  108 

2.62  x  10s 

Calculated  Maximum  Stress  without  Bending  Factor  (lb. /in.2) 

37,  560 

28,  620 

(N/m2) 

2.  59  x  108 

1.97  x  108 

Percent  Error  (%) 

15.9 

32.6 

fi 

1.20 

1.20 

y 

.328 

1.0 

Calculated  Maximum  Stress  with  Bending  Factor  (lb. /in.2) 

39,  700 

34,  300 

(N/m2) 

2.  74  x  108 

2.36  x  108 

Percent  Error  (% ) 

9.6 

10.7 

*  For  this  ring  an  effective  moment  of  inertia  and  effective  area  were  calculated  using 
the  approach  given  in  reference  9  . 


where  y  is  the  ratio  of  the  bending  stress  to  the  maximum  total  .stress .  Table  II 
shows  the  calculated  maximum  stress  for  these  cylinders  both  without  and  with  the 
bending  factor /3  .  For  the  latter,  /3  has  been  taken  as  1. 20  which  is  the,  ratio  of  the 
actual  failure  stress  of  cylinder  No.  7  to  the  actual  failure  stress  of  cylinder  No.  1. 
These  cylinders  are  identical;  one  was  tested  in  bending  and  the  other  in  compression. 

Additional  Remarks 


The  average  number  of  rings  ( a)  per  longitudinal  half  wavelength  for  the 
cylinders  given  in  Tables  I  and  n  is  relatively  low,  a  varying  from  1.  7  for 
cylinders  No.  2  and  No.  3  to  2.  6  for  cylinder  No.  6.  Van  der  Neut  [  6  ]  performed 
a  study  to  determine  what  error  is  produced  by  using  a  "smeared"  ring  approach 
when  a  is  low.  He  states  that  for  stiffened  cylinders  the  error  is  on  the  order  of 
4  percent  for  a  =  2.  0  and  6  percent  for  a  =  1.6,  the  exact  error  being  dependent 
upon  the  cylinder  stiffness  properties. 


|A11  the  cylinders  given  in  Tables  I  and  II  had  some  end  fixity.  Deflection 
measurements  indicate  that  for  cylinders  No.  4,  5,  and  6  the  amount  of  end 
fixity  was  negligible.  Unfortunately  there  were  not  sufficient  deflection  measure¬ 
ments  on  the  other  cylinders  to  determine  the  amount  of  end  fixity,  but  it  is 
believed  to  be  small  for  these  cylinders  also^j  -  p  ^  ^  'fsf'A  3*  $ 


DISCUSSION 


The  method  for  predicting  general  instability  developed  in  this  report  con¬ 
siders  the  eccentricity  of  the  ring  with  respect  to  the  skin  centerline.  It  can  be 
shown  that  this  factor  has  a  large  effect  on  the  general  instability  failure  load.  In 
fact,  moving  the  rings  from  the  inside  to  the  outside  of  the  cylinder  can  sometimes 
change  the  general  instability  load  100  percent  or  more.  As  an  example,  the  general 
instability  failure  load  of  cylinder  No.  4  is  654  kips  (2.  9  x  106  N)  for  inside  rings 
and  1254  kips  (5.  6  x  106  N)  if  the  same  rings  are  on  the  outside  of  the  corrugation. 
This  same  type  of  effect  is  present  in  cylinders  with  inside  or  outside  longitudinal 
stiffeners  as  is  shown  by  the  test  data  in  reference  4. 

As  the  test  results  show,  good  agreement  is  obtained  for  the  cylinders 
loaded  in  bending  and  for  those  loaded  simultaneously  in  bending  and  axial  com¬ 
pression,  if  a  bending  factor  of  1.  20  is  applied  to  the  bending  portion  of  the  load. 
Cylinders  loaded  in  bending  carry  a  greater  maximum  stress  because  the  portion 
of  the  cylinder  that  is  highly  loaded  is  stabilized  by  the  remainder  of  the  cylinder. 

At  the  present,  though,  sufficient  information  is  not  available  as  to  exactly  what 
bending  factor  should  be  used  for  each  particular  cylinder,  this  factor  being  a 
function  of  the  cylinder  stiffness  properties.  Based  on  currently  available  infor¬ 
mation,  it  is  not  recommended  that  a  bending  factor  be  used  for  design. 
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The  use  of  linear  small  deflection  theory  has  been  justified  because  the 
relatively  high  bending  stiffness  of  a  corrugated  shell  makes  the  cylinder  less 
susceptible  to  initial  imperfections  and  other  monocoque  effects.  A  corrugated 
shell  for  which  this  assumption  does  not  hold  true  is  not  envisioned  as  being 
practical,  but  still  it  would  be  possible  to  make  a  cylinder  with  such  a  small 
corrugation  depth  that  it  would  actually  approach  monocoque  in  properties.  For 
a  design  of  this  nature  a  reduction  factor  will  have  to  be  applied  to  the  failure 
load  given  by  this  method.  An  acceptable  technique  for  determining  the  reduction 
factor  is  given  by  Almroth  [  8]  . 

For  all  the  cylinders  tested  the  ring  extensional  and  in-plane  bending  stiff¬ 
nesses  were  the  only  ring  properties  affecting  the  failure  load.  It  is  believed 
by  the  authors  that  for  most  practical  applications  the  out-of-plane  bending- 
stiffness  and  the  torsional  stiffness  of  the  ring  are  of  secondary  importance. 

In  computing  the  stiffness  properties  of  the  ring,  care  should  be  taken  by 
the  analyst  that  the  effective  rather  than  the  apparent  stiffness  properties  are 
used.  This  is  especially  true  for  the  out-of-plane  properties  such  as  the  lateral 
bending  and  torsional  stiffness,  but  may  also  be  important  for  the  in-plane  proper¬ 
ties  of  the  ring;  e.  g.  a  channel  section  having  wide  and  comparatively  thin  flanges 
may  not  be  fully  effective  in  bending,  as  is  discussed  in  reference  9. 

It  should  also  be  mentioned  that  the  method  given  in  this  report  only  de¬ 
termines  the  general  instability  failure  load.  It  does  not  check  for  panel  buckling 
or  local  crippling  failures. 


CONCLUDING  REMARKS 


Linear  small  deflection  theory  has  been  used  to  develop  a  method  to  de¬ 
termine  the  general  instability  load  of  a  ring-stiffened  corrugated  cylinder  under 
axial  compression.  The  general  instability  failure  loads  of  seven  corrugated  cylin¬ 
ders  have  been  compared  with  loads  calculated  by  using  this  method.  Agreement 
between  the  calculated  and  actual  failure  loads  were  quite  close  for  the  five 
cylinders  loaded  in  pure  compression.  When  a  bending  factor  of  i.  20  was  used, 
good  agreement  was  also  obtained  for  one  cylinder  loaded  in  pure  bending 
and  for  one  cylinder  having  a  combination  of  axial  and  bending  load.  Since 
calculations  must  be  made  for  many  different  mode  shapes  before  the  minimum 
buckling  load  can  be  determined,  a  computer  program  was  developed. 
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FIGURE  1.  FORCES  AND  MOMENTS  ACTING  ON  SKIN  ELEMENT 
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Ring  Radial  Deflection  in.  (cm) 


FIGURE  6.  CIRCUMFERENTIAL  MODE  SHAPE  OF  CYLINDER  TESTED  IN  AXIAL 
COMPRESSION 
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Failure  Load 


APPENDIX 


The  computer  program  is  written  in  Fortran  II.  The  program  notation  is 
given  on  page  24  and  the  program  listing  on  pages  25  through  27.  The  sequence 
of  steps  for  program  compilation  and  operation  is  given  in  reference  10. 

The  program  input  data  and  corresponding  column  locations  on  the  input 
cards  are  in  the  following  form. 


1  -  19 

20  -  69 

70  -  80 

Card  1 

- 

Case  Title 

- 

1-16  . 

17-32 

33-48 

49-64 

65-80 

Card  2 

E 

ER 

G 

GR 

R 

Card  3 

QL 

TX 

TS 

QIX 

C 

Card  4 

QLR 

AR 

QIYR 

QIZR 

QJR 

1-4 

5-8 

9-12 

13-16 

17-80 

Card  5 

Ml 

MM 

Ni 

NN 

- 

where  Ml,  N1  and  MM,  NN  are  the  first  and  last  values  of  m  and  n  respectively, 
and 


TX  =  t  (D.F.) 

TS  =  t/D.  F. 

where  D.  F.  ,  the  development  factor,  is  defined  as  the  ratio  of  the  average  skin 
area  per  inch  of  circumference  to  the  skin  thickness. 
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An  example  of  the  program  output  is  shown  on  pages  28  and  29.  The 
axisymmetric  mode  shapes  (n  =  0)  have  been  included  for  completeness. 

In  using  the  computer  program  a  sufficient  range  of  mode  shapes  must 
be  considered  so  that  the  minimum  biickling  load  is  definitely  obtained.  This  is 
mentioned  because  it  is  possible  to  have  more  than  one  apparent  minimum  buckling 
load  as  is  shown  in  Figure  10. 

COMPUTER  PROGRAM  NOTATION 


TEXT  NOTATION 

PROGRAM  NOTATION 

a.. 

AIJ 

ij 

c 

C 

m,  n 

QM,  QN 

m 

QMBAR 

n 

QNBAR 

q 

Q3 

TX,  TS 

^r 

AR 

Dx>  Dyr’  Dzr 

DX,  DYR,  DZR 

E,  Er 

E,  ER 

E 

EBARX 

Er 

EBARR 

G,  Gr 

G,  GR 

G 

GBAR 

Ix>  lyr’  *zr 

QIX,  QIYR,  QIZR 

Jr 

QJR 

Kr 

QKR 

L,  Lr 

QL,  QLR 

R,  Rr 

R,  RR 
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C  THE  GENERAL  INSTABILITY  OF  RING  STIFFENED  CORRUGATED 

C  CYLINDERS  UNDER  AXIAL  COMPRESSION 

1  READ  4 

READ  2,  E»  ER»  G»  GR»  R»  QL  »  TX»  TS*  Q1X» 

1  C »  QLR »  AR*  Q I Y  R  »  GIZR»  Q  J  R  »  Ml*  MM*  Nl»  NN 
PUNCH' 3 
PUNCH  4 
PUNCH  5 
PUNCH  12 
PUNCH  6 

PUNCH  7*  E »  ER »  G*  GR*  R 
PUNCH  8 

PUNCH  9*  OL  »  TX  »  TS*  C I X  *  C 
PUNCH  10 

PUNCH  9*  QLR,  AR,  QIYR,  QIZR*  Q JR 
PUNCH  13 

PUNCH  14,  Ml,  MM,  Nl,  NN 

PUNCH  15 

PUNCH  16 

PUNCH  17 

PUNCH  18 

EBARX  =  E*T  X 

GBAR  =  G*TS 

DX  =  E*QIX 

EBARR  =  ER*AR/QLR 

DYR  =  ER*Q I YR/QLR 

DZR  =  ER*QIZR/QLR 

QKR  =  GR*GJR/QLK 

RR  =  R  +  C 

RRR  =  R/RR 

Y  =  0 

20  DO  51  N  =  Nl,  NN 
DO  51  M  =  Ml,  MM 
OM  =  M 
ON  =  N 

QMBAR  =  QM*3 .14159/QL 

QNBAR  =  QN/R 

OR  =  R*C*QNBAR**2  +1. 

ORR  =  -R*C*QNBAR**2  +  l./RRR 

All  =  -EBARX*GMBAR**2  -  G8AK*GNBAK**2  -  ( DZK*QNBAR**2  +  Qi<R/R'**2)* 
1  RRR**3*QNBAR**2 

A12  =  G  B  A  R  *  0  M  B  A  R  *  0  N  BAR 

A 1 3  =  —  ( DZR*GKk  +  GKR ) * ( KRR*QNBAR ) **2*UMdAR/ RR 

A22  =  -  GBAR*Gi>i8AR**2  -  (EBARR/RRR  +  DYR*KRR/ R**2  )  *GNBAR**2 

A23  =  -  ( EBARR*QNBAR*Gk/R  +  DYR*QNBAR**3/RR ) 

A 3 3  =  -  DX*QKE5AR**4  -  ( EBARR*GR**2/R**2  +  DYR*QN8AR**4  +  RRR**2* 

1  <DZR*QKR**2/R**2  +  QKR*ONBAR**2  >  *GM3AR**2  )  *RKR 

DD  =  QMBAR** 6 

CC  =  QMBAR**4*(AH  +  A22  +  A33) 

BB  =  OM  B  A  R  *  *  2  * ( A 1 1  *  A  2  2  +  A11*A33  +  A22*A33  -  A12**2  - 
1  A 1 3**2  -  A23**2 ) 

AA  =  A 1 1  *  A  2  2  *  A  3  3  +  2 . *A 1 2* A  1 3* A2 3  -  A23**2*A11  -  A13**2*A22 
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1  -  A12**2*A33 

XXX  =  0. 

01  =  0. 

El  =  DD*Q1**3  +  CC  *01**2  +  BB*01  +  AA 
DQ  =  . 1*AA/ ( ( A12**2  -  A1 1*A22 ) * ( QMBAR**2 ) ) 

21  Q2  =  Q1  +  00 

22  E 2  =  DD*Q2**3  +  CC*G2**2  +  BB*02  +  AA 
IF  (E1*E2)  20*28,26 

26  Q 1  =  Q2 
El  =  E  2 
60  TO  21 

28  IF  (XXX)  29,29,30 

29  DO  =  .01*02 
XXX  =  1. 

GO  TO  21 

30  Q3  =  01 

RAD  =  -  3.*(DD*G3>**2  -  2.*DD*CC*Q3  +  CC**2  -  4«*DD*3B 
IF  (RAD)  50,  31,  31 

31  04  =  (-(CC  +  DD*Q3)  +  SORT F ( R AD ) ) / ( 2 . *DD ) 

Q5  =  (-(CC  +  DD*Q3)  -  SQRTF ( RAD ) ) / ( 2 . *DD ) 

IF  (04)  34,  34,  32 

32  IF  ( Q3  -  04)  34,  34,  33 

33  03  =  Q4 

34  IF  (05)  50,  50,  35 

35  IF  (03  -  05)  50,  50,  36 

36  03  =  Q5 

50  P  =  6.28*R*Q3 
STRES  =  Q3/TX 

IF  ( Y )  55,  55,  52 
55  STREX  =  STRES 
MX  =  M 
NX  =  N 
GO  .TO  54 

52  IF  (STRES  -  STREX)  53,  51,  51 

53  STREX  =  STRES 
MX  =.  M 

NX  =  N 
GO  TO  51 

54  Y  =  1. 

51  PUNCH  19,  03,  P,  STRES,  M*  N 
PUNCH  23 

PUNCH  24 

PUNCH  11,  STREX,  MX,  NX 
GO  TO  1 

2  FORMAT  (4F16.0,  F16.4/  5F16.4/  5F16.4/  414) 

3  FORMAT  (//19X41HTHE  GENERAL  INSTABILITY  OF  RING  STIFFENED/ 
119X44H CORRUGATED  CYLINDERS  UNDER  AXIAL  COMPRESSION//) 

4  FORMAT  (19X49H  ) 

5  FORMAT  ( ///35X10HI NPUT  DATA) 

6  FORMAT  (15X1HE,  14X2HLR,  15X1HG,  14X2HGR,  15X1HR) 

7  FORMAT  (4F16.0,  F16.4//) 

8  FORMAT  ( 14X2HGL ,  14X2HTX*  14X2HTS »  13X3HQIX,  15X1HC) 


9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

23 

24 
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FORMAT  (5F16.6//) 

FORMAT  ( 1 3X3HQLR »  14X2HAR*  1 2  X4HQ I YR  »  12X4HQIZR*  13X3HQJR) 

FORMAT  ( 4XF 1 6 ♦ 4 *  2(7X13)) 

FORMAT  ( 35X10H///////////// ) 

FORMAT  ( 1 4X2HM1 *  14X2HMM*  14X2HN1*  14X2HNN) 

FORMAT  (4(12X14)///) 

FORMAT  ( 35X11HOUTPUT  DATA) 

FORMAT  (35X11H//////////////) 

FORMAT  (15X5HAXIAL) 

FORMAT  (7X13HLOAD  PER  INCH*  10X10HAXIAL  LOAD  »  8X12HAXIAL  STRESS* 
1  9X 1 HM  »  9X1HN/ ) 

FORMAT  ( 3 (4XF16.4) *  2(7X13)) 

FORMAT  (///3X17HMIN.  AXIAL  STRESS) 

FORMAT  (2X18HIN  THE  ABOVE  RANGE*  9X1 HM  *  9X1 HN ) 

FORMAT  (F16.4) 

END 
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the  general  instability  of  ring  stiffened 

CORRUGATED  CYLINDER'S  UNDER  AXIAL  COMPRESSION 
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